NEW RESULTS ON QUANTUM CHAOS IN ATOMIC NUCLEI 



V.R. MANFREDI 

Dipartimento di Fisica "G. Galilei" 
Umversita di Padova, Via Marzolo 8, 1-35131 
Istituto Nazionale di Fisica Nucleare 
Sezione di Padova, Via Marzolo 8, 1-35131 Padova, Italy 

L. SALASNICH 

Dipartimento di Matematica Pura ed Applicata 
Universita di Padova, Via Belzoni 7, 1-35131 Padova, Italy 

Istituto Nazionale di Fisica Nucleare 
Sezione di Padova, Via Marzolo 8, 1-35131 Padova, Italy 

In atomic nuclei, ordered and chaotic states generally coexist. In this paper the 
transition from ordered to chaotic states will be discussed in the framework of 
roto— vibrational and shell models. In particular for 160 Gd, in the roto— vibrational 
model, the Poincare sections clearly show the transition from order to chaos for 
different values of rotational frequency. Furthermore, the spectral statistics of 
low-lying states of several fp shell nuclei are studied with realistic shell— model 
calculations. 



1 Introduction 

In atomic nuclei, just as in other many-body systems, ordered and chaotic 
states ^ generally coexist. In fact, in zero-order approximation the relevant 
elementary excitations (such as rotation, vibration and single-particle) may 
be regarded as independent modes. Then we consider the interaction between 
these elementary modes 1 ). At one end of the chain of complexity we have a 
single mode which can be considered " regular" , whereas at the other end there 
are the so-called "stochastic" or "chaotic states". 

In order to discuss the coexistence in atomic nuclei of ordered and chaotic 
states, many models have been used 2 - 1 ' 3 ) but, for reasons of space, we shall 
mention only the roto-vibrational and shell models. 



a Here we use the term chaos for a quantal system in a very restricted sense: " Quantum 
Chaology is the study of semiclassical behaviour characteristics of systems whose classical 
motion exhibits chaos". If the semiclassical limit is not known we use spectral statistics 4 ) 
(e.g. P(s) and As(L)) to distinguish between ordered and chaotic states. 
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2 The Roto Vibrational Model 

This model has been amply described in Ref. 5, 6 and 7 and we limit ourselves 
to reporting only a few basic formulae. The hamiltonian is 

H = H vlb + H rot , (1) 



where 



H vib = -B(al + 2a 2 2 ) + V(a ,a 2 ), (2) 



1 3 

Hrot = 2 ^2 UJ k J k( a o,a2), (3) 



fe=i 



with 



V(a ,a 2 ) = -C 2 (a 2 + 2a 2 ) + x —C 3 a (6a 2 2 -a 2 ) + -C 4 (a 2 + 2a 2 2 ) 2 + V . (4) 
z V oo 5 

The shape of the nuclear potential V(ao, a 2 ) is function of C 2 and \ = C 3 /(C 2 C4) 7 ' 1 . 
The parameters ao and a 2 are connected to the deformation (i and asymmetry 
7 by the relations 

a = /3cos7, a 2 = ^y=sin7. (5) 
In terms of the new variables the components of the moment of inertia are 8 ) 

J k =4Bf3 2 S m 2 (j-—k). (6) 
If the nucleus has an axially symmetric equilibrium configuration and 

u 1 =u 2 = -j=, w 3 = 0, (7) 
the hamiltonian (1) can be written as 

H = \B{a 2 + 2a\) + V(a , a 2 ) + ^Btu 2 (3a 2 + 2a 2 2 ), (8) 

where V(ao, a 2 ) is given by (4) and Vo is chosen to have the minimum of the 
potential equal to zero. 



Figure 1: The Poincare sections for 160 Gd at the energy 5.5 MeV and for different val- 
ues of rotational frequency; from the top: hw = MeV, hu> = 0.5 MeV, hui = 1 MeV. 
Adapted from Rcf. 9. 
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3 Numerical study of the order chaos transition 

A very useful tool for the study of the global instability is provided by the 
Poincare sections. The classical trajectories have been calculated by a fourth 
order Runge-Kutta method 9 ' . 

The Hamilton equations are as follows 

d = Bp a , 

d 2 = 2Bp 2 , 

I 2 4 
Po = -C 2 a - 2d — C 3 (3aj - 3a„) - -C 4 a (al + 2a 2 2 ) - 3Buj 2 a a , (9) 

/ 2 g 

p 2 = -2C 2 a 2 - 12 W — C 3 a a 2 - -C 4 a 2 (al + 2o|) - 2Buj 2 a 2 , 
where po and p 2 are the conjugate momenta 

po = Bd Q , p 2 = 2Bd 2 . (10) 

Figure 1 shows the Poincare sections for 160 Gd at the energy 5.5 MeV 
and for different values of rotational frequency. The figure clearly shows a 
chaos-order transition as the frequency w increases. In Figure 2, for 166 Er, 
the Poincare sections are shown for different values of the energy and rotational 
frequency lo — 0. As can be seen, there is a chaos-order transition, but not so 
sharp as in the previous case. 

Fluctuation properties of quantal systems with underlying classical chaotic 
behaviour and time-reversal symmetry are in agreement with the predictions 
of the Gaussian Orthogonal Ensemble (GOE) , and quantum analogs of classi- 
cally integrable systems display the characteristics of Poisson statistics 10 '. In 
general, various statistics may be used to show the local correlations of the 
energy levels; we shall discuss P(s) and A 3 (L) only. P(s) measures the prob- 
ability that two neighbouring eigenvalues are a distance "s" apart. For GOE 
we have the Wigner distribution 

P(s) = ^seM~s% (11) 

which gives level repulsion. A 3 (L) is defined for a fixed interval (— L/2,L/2) 
as the least-square deviation of the staircase function N(E) from the best 
straight line fitting it 



1 r L/2 

A 3 (L) = -mm [N(E) - AE - B] 2 dE, (12) 

L A.B J _ L / 2 
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Figure 2: The Poincare sections for 166 Er at the rotational frequency ui = and for different 
values of the energy: (a) E = 1 MeV, (b) E = 6 MeV, (c) £ = 9 MeV, (d) E = 12 MeV. 
Adapted from Ref. 9. 
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Figure 3: Spectral statistics P(s) and Aa(L) for 166 Er at the rotational frequency oj = 
for different energy regions: 2 < E < 6 MeV (below) and for 13 < E < 17 MeV 
(above). The solid line is the GOE statistic curve and the dashed line is the Poisson one. 
Adapted from Ref. 9. 



where N(E) is the number of levels between E and zero for positive energy, 
between —E and zero for negative energy. A^(L) provides a measure of the 
degree of rigidity of the spectrum: for a given interval L, the smaller A^,(L) 
is, the stronger is the rigidity, signifying the long-range correlations between 
levels. For this statistics in the GOE ensemble 



In Figure 3 the spectral statistics P(s) and A 3 are plotted for 166 Er. These 
statistics confirm the classical results: for energies above the saddle energy (<~ 4 
MeV) there is prevalently chaotic behaviour; for higher energies there is mixed 
behaviour. 




(13) 
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4 Shell Model Calculations 



In this section we discuss the statistical analysis of the shell-model energy 
levels in the A = 46-50 region. Exact calculations are performed in the 
(/7/2>P3/2)/5/2)Pi/2) shell-model space, assuming 40 Ca as an inert core 11 ). Di- 
agonalizations are performed in the m-scheme using a fast implementation of 
the Lanczos algorithm with the code ANTOINE. For a fixed number of valence 
protons and neutrons we calculate the energy spectrum for projected total an- 
gular momentum J and total isospin T. The interaction we use is a minimally 
modified Kuo-Brown realistic force with monopole improvements. 

We calculate the T = T z states from J = to J = 9 for all the combina- 
tions of 6 active nucleons, i.e. 46 V, 46 Ti, 46 Sc, and 46 Ca, and also for 48 Ca and 
50 Ca. 

Since we are looking for deviations from chaotic features, we are mainly 
interested in the low-lying levels, up to a few MeV above the JT yrast line. 
Let us consider the energy levels up to 4, 5 and 6 MeV above the yrast line, 
and calculate the fluctuations around the average spacing between neighbor- 
ing levels. In this range of energies, the level spectrum can be mapped into 
unfolded levels with quasi-uniform level density by using the constant tem- 
perature formula. To guarantee that the results up to different energies are 
unaffected by the unfolding procedure, the unfolding is performed using always 
the whole set of levels up to 6 MeV, for each JT set in the nucleus. 

The mean level density can be assumed to be of the form 

p(E) = ±e X p[(E-E )/T}, (14) 

where T and Eq are constants. For fitting purposes it is better to use not p(E) 
but its integral N(E). We write 

N(E) = / p(E')dE' + N = exp \{E - E )/T] - exp [-E /T] + N . (15) 
Jo 

The constant No represents the number of levels with energies less than zero. 
We consider this function as an empirical function to fit the data and let Nq 
take non-zero values. The parameters T, Eq and No that best fit N{E) are 
obtained by minimizing the function 

G(T,E ,N )= / [N(E) - N(E)] 2 dE, (16) 

where N{E) is the number of levels with energies less than or equal to E. The 
energies E min and E max are taken as the first and last energies of the level 
sequence. 
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Figure 4: P(s) distribution for low— lying levels of fp shell nuclei with < J < 9: (a) 46 V, 
46 Ti and 46 5c; (b) 46 Ca, 48 Ca and 50 Ca. The dotted, dashed and solid curves stand for 
GOE, Poisson and Brody distributions, respectively. Adapted from Ref. 11. 



As previously discussed, the spectral statistic P(s) may be used to study 
the local fluctuations of the energy levels. P(s) is the distribution of nearest- 
neighbour spacings Si = E i+ i — Ei of the unfolded levels. 

For quantum systems whose classical analogs are integrable, P(s) is ex- 
pected to follow the Poisson limit, i.e. P(s) = exp (— s). On the other hand, 
quantal analogs of chaotic systems exhibit the spectral properties of GOE with 
P( S ) = (7r/2) S exp(~f S 2 ). 

To quantify the chaoticity of P(s) in terms of a parameter, it can be 
compared to the Brody distribution, 

P(s,uj) =a(w + l)s w cxp(-as w+1 ), (17) 

with 

a = (T[^}r+\ (18) 

LO -f- 1 

This distribution interpolates between the Poisson distribution (ui = 0) of 
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intcgrable systems and the GOE distribution (ui = 1). The parameter u) can 
be used as a simple quantitative measure of the degree of chaoticity. 

The number of J = 0-9 spacings below 4, 5 and 6 MeV range from 42, 66 
and 105 in 46 Ca, to 86, 149 and 231 in 46 Ti, respectively. 

To obtain a better estimate of the Brody parameter, we can combine spac- 
ings of different nuclei. The number of level spacings is now sufficiently large 
to yield meaningful statistics and we see that Ca isotopes are not very chaotic 
at low energy, in contrast to other nuclei in the same region (Figure 4). 

Why are Ca isotopes less chaotic than their neighbors? We observe that 
the two-body matrix elements of the proton-neutron interaction are, on aver- 
age, larger than those of the proton-proton and neutron-neutron interactions. 
Consequently the single-particle mean-field motion in nuclei with both pro- 
tons and neutrons in the valence orbits suffers more disturbance and is thus 
more chaotic. 

5 Conclusions 

We have shown that in the roto-vibrational model of atomic nuclei, an order- 
chaos-order transition occurs as a function of the energy. Concerning the shell 
model calculations, the main conclusion of this paper is that for Ca isotopes we 
find significant deviations from the predictions of the random-matrix theory. 
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